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6o '. Abstract 

The interaction theory between charged scalar field and electromagnetic field on de Sitter back- 
ground is studied in the hypothesis of minimally coupling. Then the reduction rules for scalar field 
"^ ■ are established and further the theory of scattering for scalar field is constructed using perturba- 

fvi . tion theory. With these elements the scattering of a charged scalar field on Coulomb potential 

is studied. One obtains that in contrast with what happens on Minkowski case the modulus of 
momentum is not conserved in the scattering process on de Sitter space. After that we study the 
problem of calculating the cross section, pointing out a series of difficulties that arises. 
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I. INTRODUCTION 

Quantum field theory on de Sitter space-time was intensively studied. The reason for this 
special attention is due to the fact that de Sitter space is the unique maximally symmetric 
curved space-time. It enjoys the same degree of symmetry as Minkowski space (ten Killing 
vectors). Even with this facilities the presence of curvature and non-trivial global properties, 
introduce new problems related to the quantization of fields on de Sitter space. Also it turned 
out that the de Sitter metric played a central role in the inflationary cosmological scenario. 
These reasons explain the enormous amount of literature on de Sitter quantum field theory. 

The quantum scalar field is studied in problems concerning quantum effects in presence 
of gravitation and cosmology. Of a special interest in cosmology is the de Sitter expanding 
universe carrying scalar fields variously coupled with gravitation. In the quantum theory a 
special role is played by the free fields (minimally coupled with gravitation) because with 
their help we can calculate scattering amplitudes using perturbation methods. In our opinion 
the calculations of scattering processes must be the principal goal of the quantum field theory 
in the curved space. This implies a theory of scattering on curved spacetime which must be 
obtained using analogue methods as in flat space case. 

The theory of quantum scalar field on de Sitter space was intensively studied, but the 
explicit solutions was written in a few papers. In our opinion these solutions have a main 
importance, in developing the theory of interaction on de Sitter space-time. A detailed study 
of interaction between quantum fields minimally coupled with gravitation in a expanding 
universe was not done from our knowledge. We focus our attention on this subject because 
the last astronomical observations shows that the expansion of the Universe is accelerating, 
thus increasing the interest on de Sitter model. The majority of investigation related to the 
quantum field theory on curved space-times it concerns problems related to the creation of 
particles on curved space, but neglect the scattering processes. This theory of scattering 
could be of interest, since the study of interactions between fermions, bosons and scalar 
fields on a curved background will help us to understand how the well established picture 



from flat space wil 



be modified by a gravitational field. Until now detailed studies related to 



the Dirac field [2[], 3j , [13( , [14| and scalar field [9| on de Sitter space-time was done obtaining 
important results such as normalized analytical solutions and quantization. The less studied 
was the electromagnetic field, but here the situation is simple is one recalls the conformal 



invariance between de Sitter and Minkowski metrics. Thus at this moment seems that we 
have all elements for constructing a theory of interaction between quantum fields on de Sitter 
background. In this paper our attention will be focused on the interaction between scalar 
field and electromagnetic field. 

We start in the second section with a brief review of the main results concerning the 
quantum scalar field theory on de Sitter space. In section IIHI we construct the theory of 
interaction between scalar field and electromagnetic field on de Sitter space and then our 
attention will be focused on developing the reduction formalism for scalar field. Section 
HVl is dedicated to the perturbation theory which will help us to obtain the definition of 
scattering amplitudes. The scattering process that implies the scattering of a charged scalar 
particle on Coulomb potential on de Sitter space is studied in section [V] and in section ED 
we study the problems that arise when we pass from scattering amplitude to cross section. 
Our conclusions is summarized in section IVIH and in Appendixes we give the main steps 
leading to our scattering amplitude and the form of the integrals that help us to establish the 
definitive form of the cross section. These results are presented in natural units h — c — 1. 

II. PLANE WAVE SOLUTIONS FOR SCALAR FIELD 

We start with the exact solutions of the free Klein-Gordon equation on de Sitter space- 
time, which was obtained in |9(,[15|. Let us write the de Sitter line element (lj: 



where u is the expansion 



Gordon equation reads 15] 



ds 2 = dt 2 - e 2uJt dx 2 } (1) 

actor. In the chart {t, x} with Cartesian coordinates the Klein- 



(8 2 - e~ 2uJt A + 3uod t + m 2 ) <j)(x) = . (2) 

It is known that the Klein-Gordon equation (T5]) on de Sitter space can be analytically solved 

in terms of Bessel functions J9|],[10J]. From different versions of solutions which are currently 

used, we will use here the normalized solutions of positive frequencies which can be wreathed 

as [lj: _ 

1 [W e-^l 2 _^ to „m/p 



''<*>= sV^'"**^ £'"**) ^ < 3 > 



where H^z) is a Hankel function of first kind, p = \p\ and the notations are: 

ni 



k = \l^ 2 - ! » ^ 



with ?7i > 3a;/2 and the fundamental solutions o 
solutions satisfies the ortonormalizations relations 



oo 



(4) 



negative frequencies are fff(x). These 



m 



ij d 3 x(-g) 1/2 fi(x)d t f pi (x) = 
-tfdPxi-gYPffWdtfrix) = 5 3 (p-p') 
% [d 3 x(-g)^ 2 f f (x)d t fp'(x) = 0, 



(5) 
(6) 



where the integration extends 
the completeness condition 



on an arbitrary hypersurface t = const and (— g) 1 ^ 2 = e 3wi ,and 
15]: 



i / d s p f}(t,x) d t fa(t,x') = e-^5\x - x') 



(7) 



The quantization can be done considering the plane wave in momentum representation: 

0(x) = 0«(x) + 0(")(x) = j d 3 p [f f (x)a(p) + f^(x)b + (p)] , (8) 

where the particle (a, at) and antiparticle (b, &') operators satisfy the standard commutation 
relations in momentum representation: 

[a(p), a\p')} = [b(p), b\p')} = 5 3 (p- f) . (9) 

These operators can be calculated using the inversion formulas 

a(p) = ijd'xi-gfl'f^dt^x) 

b{p) = ijd'xi-g) 1 ' 2 ^) 5 ^(x), (10) 

from which the formulas for a\tf can be obtained immediately. All these operators act on 
the Fock space supposed to have an unique vacuum state |0) accomplishing 



a(p)|0) = 6(p)|0) = , (0\a\p) = (0\b j {p) = . 



(11) 



In this way the field <ft is correctly quantized according to the canonical rule [15] 



[<p(t,x),d t tf(t,x')] =ie- 3uJt 5 3 (x-x') 



(12) 



In the quantum theory of fields it is important to study the Green functions related to the 
partial commutator functions (of positive or negative frequencies) defined as: 

D (±) (x,x') = i[</> (±) (x),0 (±)t (a/)L (13) 

and the total one, D = D^ + D^~\ These function are solutions of the Klein-Gordon 
equation in both the sets of variables and obey [D^(x, x')]* = D^(x,x') . Then the 
partial commutator functions will have the expressions: 

D<> + \x,a?)=iJ<Fpf f (x)%&), 
DH(z,a/) = -ijd 3 pft(x)f p (x'). (14) 

These functions help us to introduce the principal Green functions. In general, G(x, x') = 
G(t, t', x — x ') is a Green function of the Klein-Gordon equation if satisfies [151 ]: 

(d 2 - e~ 2uJt A x + 3cu<9 t + m 2 ) G(x, x') = e^d^x - x') . (15) 

The properties of the commutator functions allow us to construct the Green function just 
as in the scalar theory on Minkowski spacetime. We assume that the retarded, Dr, and 
advanced, D A , Green functions read: 

D R (t, t',x-x') = 6(t - t')D(t, t',x-x'), (16) 

D A (t,t',x-x') = -6{t' -t)D(t,t!,x-x'), (17) 

while the Feynman propagator, 

D F (t,t',x-x') = z(O|T[0(ar)0 t (x / )] |0) 

= e{t-t')D^\t,t\x-x')-e{t' -^D^it.t'.x-x'), (18) 

is defined as a causal Green function. Thus we have now all the elements that we need for 
begin our study related to the reduction formulas for scalar field on de Sitter space-time. 

III. REDUCTION FORMALISM 

As one know the goal of the quantum field theory is to describe the dynamics of interacting 
particles which are observed in nature. For that we are interested to obtain the amplitudes 



for single particles to propagate in a curved space-time and the transition amplitudes for 
interacting particles between different initial and final states. Our main attention in this 
paper will be focused on interaction between charged massive scalar field and electromagnetic 
field on de Sitter space-time, because a theory of free fields can't give us details about 
physical quantities that can be measured. In this work we use for our calculation the same 
formalism as in [4] . [51] . [1 1| 



121 ]. and we adopt the minimal coupling corresponding to the 
classical interaction of a point charge as the prescription for introducing electrodynamic 
couplings. Also we will note the interacting fields with <f>{x) and A^x). These fields will 
satisfy one system of equations that can be obtained from an action of the free scalar field, 
free electromagnetic field and an interaction term, all minimally coupled with gravitational 
field: 

-ie[{d^)<f> - (^0)0^ + eW^ \, (19) 



^ - m 2 0V) - | V 



where the density of lagrangian of interaction between the scalar charged field and electro- 
magnetic field d = —ie^f—g^d^)^ ~~ (9 IM cj))cj)^]A fl + ^f—ge 2 '$ (frA^A^ was obtained after 
we make the substitution <9^0 — > d^cj) — ieA^ip in the lagrangian of scalar field. From action 
l fl9l) . we obtain the system of equations that describe the interactions between the charged 
scalar field and electromagnetic field on a curved background as follows: 



y=d,[^g)F^) = -ieWtpW-p 



+2e 2 (p ] <j)A h 
1 



<9 M y=ggTd v 4> - ie^^A*] + m 2 = ie{d^)A^ + e 2 <j>A^, 
' ■ dp [y/^gg^dvfi + ie^g^A^] + m 2 tf = -ie(d^)A^ + e^U^. (20) 



The problem of constructing exact solutions for the system (1201) . which contain coupled 
nonlinear equations for interacting fields proves to be as in the flat space case, a impossible 
task. It is immediate from here that for constructing solutions of the coupled equations, 
we can try to follow the same approximation method as in flat space case. In the method 
of approximation from flat space case the system of coupled equations is replaced with one 
system of integral equations that contain information about initial conditions. We follow 
here the same steps and for that we select one Green function G(x,y) corresponding to one 



initial condition, which help us to write the solution of the second equation from system 
(120]) like follows: 



0( x ) = 0(x) + ey d i y^G(x,y)t-j=d lx [^ ( j ) (y)A»(y)] 

+i[d»<j>(y)]A»(y) + e<f>(y)A fl (y)A»(y)\, (21) 

where <f)(x) are free fields. One can verify that (12"T1) is the exact solution of the second 
equation from (1201) . applying (ExG+ m2 ) to (p(x), with the observation that [ExG+ rn ' 2 ] ( t ) (x) = 
0, where we introduce the following convenable notation Ekg = 7= d^yf^g g^ u d v ). Using 
the fact that Green functions must satisfy one equation of the type ffl5l) it is not hard to 
verify that our solution (I2TI) is justified. 

Since our main scope in this paper is to study scattering problems, we want to construct 
the states which give a simple description of the physical system at the time t — > —00. 
At this time the particles in a scattering process have not yet interacted with each other 
and propagate under the influence of their self-interactions only. This implies that for 
constructing a theory of interactions, we must seek for operators that create independent 
particle states with each particle propagating with its physical mass. Let us start to construct 
these operators, using (J2T1) which offers us the possibility of constructing free fields, which 
are asymptotic equal (at t — > ±00) with solutions of system ( |20l) . As we know the retarded 
Green functions Dr(x, y) vanishes at t —* — 00 while the advanced one Da(x, y) vanishes for 
t — > 00. Thus (I2TJ) can be expressed with retarded and advanced functions as follows: 

0(z) = 4>r/a(z) +e d A y*J=gD R/A {x,y)\ -=d fl [i/=g<j>(y)A' t (y)] 



+i[d^(y)]A»(y) + e<f>(y)A^y)A»(y)j, (22) 

then they would define free fields, 4>r/a(%) that satisfy: 

lim (<f>(x) - 4> R/A (x)) = 0. (23) 

The free fields <pR and <pA have mass m and are equal at t ± 00 with exact solutions of the 
coupled equations and represent the fields before and after the interaction. It is known that 
the mass of the interacting fields may differ from the mass of free fields. For obtaining the 
mass correction we will ad to the lagrangian of interaction on term of the form —8m 2 ' 



which help us to rewrite the second equation from system (1201) : 
1 



-g 



<9 M [V^gg^d^ - ie^g^A* 1 ] + ™ 2 = ie(d„<f>)A li + e 2 (j)A^ - 5m 2 0, (24) 



where the current j = ie(d ^0) A^ + e 2 <f)A fl A fl — 5m 2 cf) + -j^d^^—gipA^) is now the source 
giving rise to the scattered waves. 

Further we ca follow the method from Minkowski case, where the free fields cj)R and 4>a 
are defined up to a normalization constant noted with ^fz. This allows us to define the 
in j out fields in order to construct the states of single particles, as follows: 

y/z<t>in/out(x) = <f)(x) -e / d A y\f=gD R/A {x,y)\ -=d IJ ,[s/ : ^g(f)(y)A fl (y)] 



+i[d^(y)]A^y) + e<Ky)AMA»(y) - Sm 2 <f>(y)j. (25) 

The in/ out free fields defined above satisfy Klein-Gordon equation, and can be written with 
the help of creation and annihilation operators and in addition satisfy conditions: 

lim (0(x) - ^/z(f)in/out{x)) = 0. (26) 

t— >=poo 

For our further calculations it will be useful to express this fields with the help of equation 
(1241) . as follows: 

y/z4>in/out{x) = 4>{x) - / d i y^/^gD R/A (x,y)[E KG (y) +m 2 }<j)(y). (27) 

From the above normalization we can define the creation and annihilations operators as in 
Section II , 



b\p)in/out = i / d 3 x(-g) 1/2 fjj(x) d t (pin/outix). (28) 



The creation and annihilation operators defined above satisfy the commutation relations Q 
and from that it follows that all the properties of free fields will be preserved. 

Before starting our calculations we make a few remarks about the scattering operator. 
Denoting the vacuum state |0) (which is suppose to be unique), the one particle/antiparticle 
states can be written: 

aKp)in/out\0) = \in/out,l(p)) 

b\p)in/ouM = \in/ out, 1(f)). (29) 



If we consider two states \in, a) and \out, (3), then the probability of transition from state a 
to state (3 is defined as the scalar product of the two states: (out, /3\in, a). These are just 
the elements of matrix for scattering operator Sp a , which are of physical interest because 
measures the transitions that could be observed experimentally. This operator assures the 
stability of the vacuum state and one particle state, and in addition transform any out field 
in the equivalent in field. 

Our remaining task is to construct and study, general matrix elements which describe 
the dynamical behavior of interacting particles. As in Minkowski theory we can construct 
n— particle in and out states as in the free scalar field theory by repeated application to the 
vacuum of a)(p) in/out and b ] (p) in/out . 

We have now all the ingredients for starting our calculations related to the reduction of 
scalar particles on de Sitter space-time. Let us consider the amplitude of one process in 
which particles from in states denoted with a, together with a scalar particle l(p) pass in 
out state, in which we have a scalar particle l(p ') and particles denoted by (3. The amplitude 
of this process can be written as: 

(out/3, l(p')\ina, l(p)) = (out(3\a out {p')\ina, l(p)) 

= (out(3\a in (p')\ina,l(p)) 

+ (out(3\[a out (p') - a in (p')]\ina, l(p)). (30) 

The first term in ( |30l) give if one use the commutation relations ([9]) the amplitude of one 
process where the scalar particle passes from state in in out state without interacting with 
other particles: 

(out(3\ai n (p')\ina, l(p)) = 5 3 (p — p ')(outf3\ina). (31) 

The second term in ( |30l) is of physical interest and we must evolve the difference: 

a ou t(p') - a in (p') = i / d 3 xe 3uJt f^,(x) d t (<j) ut(x) - 4>in(x)). (32) 

For evolving ( !3~2l) we use equation ( |27|) which give us the in/ out fields as function of advanced 
and retarded Green functions: 



<Pout{x) - <f>i n (x) = —= I d y^/=gD(x,y)[E KG (y) + m ]<f>(y), (33) 

where we noted D(x,y) = D R (x,y) — D A (x,y). If we introduce now (1331) in (1321) and in 



addition we use the following relation obtained after a simple calculation: 

i J tfxe^fyix) 5 D(x,y) = ift,(y), (34) 

we finally obtain: 

(out/3, l(p')\ina, l(p)) = 5 3 (p — p')(out/3\ina) 

+ 4= [^f*(y)[E KG (y)+m 2 } 



x(out(3\(j)(y)\ina,l(p))d A y. (35) 

Let us apply this method for reduction of particles from in state. As above we start as 
in (130]) . but this time we apply out method to the scalar particle from in state: 



(out/3, l(p')\ina, 1(f)) = (out/3, l(p')\a\ n (p)\ina) 

= (outf3,l(p')\a ] out (p)\ina) 

+ (out(3, l(p')\[a\ n (p) - al ut (p)}\tna). (36) 

The first term in (136]) give us if we use commutation relations (Q the amplitude of one 
process where the scalar particle don't interact with other particles, and the difference from 
second term can be evolved using the method presented above: 

4niP) - oL(p) = -i fd 3 xe^f p (x) 5 (<t>l(x) - ^ out (x)), (37) 



which can be evolved using the hermitic conjugate of (p7|) and (|34|) . Finally the amplitude 
can be wreathed as: 



(out/3, l(p')\ina, l(p)) = 5 (p — p')(out/3\ina) 



i 



+^ / y/=g{autp,\ip ')W{y)\ina,l(p)) 



x[E KG (y)+m 2 }f 1 ,(y)d 4 y. (38) 

The reduction rules for antiparticles from out/in states can be obtained in the same 
manner and we give just the final expressions: 

(out(3,l(p')\ina,l(p)) = 5 3 (p — p')(outp\ina) 

x[E KG (y) + m 2 ]f f ,(y)d 4 y, (39) 

10 



{out/3, l(p')\ina,l(p)) = 5 3 (p — p'){out(3\ina) 



i 



2i 



f=-gfZ{y)[E KG {y)+m 
x {out/3, l(p')\(p(y)\ina)d 4 y. (40) 

For reduction of the second particle we use the results and the methods from reduction 
of first particle. We suppose that we already done the reduction of the first particle from 
out state and we obtain (j35p . Now if we suppose that in the particles denoted by (3 we have 
another scalar particle then (3 = (3' + l(p"). The matrix element in which we are interested 
appears in fl35l) : 

{out(3\(j){y)\ina,l{p)) = {out/3', l(p")\(p(y)\ina, l(p)) 

= {out/3' \[a out {p")(/){y) - (f>(y)a in {p"))\ina, l(p)) 

+{out/3'\(j ) (y)a in (p")\ina, 1(f)). (41) 

The last terms is the form {out/3'\</)(y)\ina)5 3 (p" — p) , and corresponds to a process where 
the scalar particle don't interact with other particles and for that reason this term is not 
interesting for us. The first term in (}4"I]) give the amplitude that interest us, and the above 
difference can be wreathed: 

a out (p")(f>(y) - (/>(y)a in (p ") = i / e 3ujt f^,,(x) d t [(p out (x)(f)(y) 

-0(y)<Mx)] d 3 %- (42) 

It is not difficult to see that equation (j27|) allows us to replace the in/out fields with expres- 
sions that depends on retarded and advanced Green functions, and using the explicit form of 
the Klein-Gordon operator in the chart {t, x}, written in Eq.Q Ekg — &t ~ e~ 2uJt A + 3udt 
, we can express our amplitude in terms of time-ordered product of field operators: 

{outp'Haoutip")^) - <j){y)a in {p")]\ina,l{p)) = -y= / V=y~ft„(z)[E KG {z) +m 2 ] 

x {out(3'\T[(j)(z)(j){y)\ \ina, l{p ))rf 4 z(43) 

From the above calculations one sees that the reduction method for the second particle 
is done using the same method as for reduction of the first particle. In the matrix element 
that interest us (14*31 . appears two field operators that are multiplied in chronological order. 
Repeating the calculation for other particles from in and out states one observe that the 
reduction calculus is the same. 

11 



Also we can obtain a generalization of the above formulas. Supposing that we have n—out 
and m — in , scalar particles, after applying the above reduction formalism to remove all 
particles from the in /out states, we arrive at the vacuum expectation value of a chronological 
product of field operators: 

(• \ m+n n „ m „ 

~r) Tl / d4x iV-9(xi) n / d A y jS J-g(y j ) 

xf^Xi^EKoixi) + m 2 ](0\T[^x i )... ( p(x n ) 
4>\y 3 )..4\y m )m[E KG { yj ) +m 2 ]f pl ( yj ). (44) 

Let us resume the above results related to reduction method. One can show that all 
particles, will be replaced by formulas that was obtain in the reduction of one particle. 
Also when more particles are reduced, in matrix element appears the time ordered products 
of corresponding field operators. Every particle will be replaced after we apply the above 
reduction method with expressions which depend on field operator <f>(x). Through reduction 
method the amplitudes was expressed as function of fundamental solutions of the free Klein- 
Gordon equation on de Sitter space-time, and as function of vacuum expectation value of 
time ordered product of fields. The vacuum expectation value of time ordered product of 
fields, together with the normalization constants {yfz) will define the Green functions of the 
interacting fields. As we know from flat space case, one can associate one Green function to 
any process of interaction that must be calculated using perturbation theory. 

In the end of this section we write the reduction rules for particles and antiparticles, 
denoting one scalar particle by l(p) and one antiparticle by lip), after reduction the in/ out 
particles will be replaced by: 



l{p') out - — J {-g)V*%,{x)[E KG (x) + m 2 ]<P(x)d 4 x, 

T(p") out _> -' [(-gy/^( x )[E KG (x)+m 2 }f f ,(x)d 4 x, (45) 



lip) in - -^=ji-g) 1/2 ^ix)[E KG ix)+m 2 ]f f ix)d 4 x, 

T(£) in -► --j= ji-g) l/2 f;ix)[E KG ix)+m 2 } ( pix)d i x. (46) 



12 



IV. PERTURBATION THEORY 

Our task here is to expand transition amplitudes and matrix elements in a power series 
in the interaction strength and obtain rules for computing the terms in the expansion. As 
we point out in the previous section the Green functions of the interacting fields can't be 
calculated exact and for that reason we must use perturbation methods. The form of the 
amplitudes obtained from reduction formalism allows us to use perturbation calculus as in 
Minkowski theory. 

For calculating the Green functions it will be helpful to use the methods from flat space 
case and write them as functions of free fields, because we know their form and properties. 
Then one Green function can be wreathed in generally as follows: 

G( yi ,y 2 ,...,y n ) = —^(0\T[^yMy 2 )..4(y n ),S]\0), (47) 

where S is a unitary operator , and have a closer form with the same operator from Minkowski 
theory and can be correlated with the scattering operator. Then entire perturbation calculus 
is based on development of operator S: 

S = Te^fW* = 1 + J2 ^J- f T[£>i{xi)~Mxn)\d A x x ...d A x n , (48) 

n=l ' ^ 

where the density lagrangian of interaction reads (we neglect the second term in the density of 
lagrangian of interaction): Cj(x) = —ie^—g : [(9^0^)0 — (d^^cfrjA^ :. Each term from (148]) 
corresponds to a rank from perturbation theory. Replacing ( 1481) in the expression of Green 
function (j47p one obtain perturbation series which allows one to calculate the amplitude in 
any order and we can write: 



(0\T[<j ) (y 1 ) ( j ) (y 2 )... ( j ) (y n ),Ci(xi)...C I (x n )]\0)d'x 1 ...d' i x n . (49) 

n!(0|ib|0) 

The evaluation of the integrant from (j4"9"l) is the same as in Minkowski case, we have a 
cinematic part which is obtained from reduction formalism and a dynamic one represented 
by operator S. The remaining task is to make the T contractions between cinematic and 
dynamic parts. 

With this method one can obtain the definition of scattering amplitudes in any order 
of perturbation theory. The definition for scattering amplitudes will help us to calculate 
scattering processes that implies scalar particles, on de Sitter expanding universe. Let 

13 



us obtain the scattering amplitude in the first order of perturbation theory. We consider 
one amplitude of the form (out,l(p')\in,l(p)) and replace the in/out particles using the 
reduction formalism developed in section IIII1 

(out,i(p')\in,i(p)} = s\p-p')-- z J f \Tg¥)V^g^)f^(y)[EK G (y) + rn 2 } 

x{0\T[cj>(y)tf(z)}\0)[E KG (z)+m 2 }f f (z)d 4 yd 4 z. (50) 

Using (I49p and neglecting the first term in (1501 . one can obtain the scattering amplitude in 
first order of perturbation theory as: 



x(0|T[: <f>{y)<ft(z) :: [(d^ (x))<f>(x) - (d^(x))^(x)]A^x) : 



x[E KG (z) +7nr]fp{z)drycrzdrx. (51) 

After performing the T contractions between cinematic and dynamic part using the method 
from Minkowski theory, and use the fact that the operator (Ekg + m2 ) acts on propagators 
as in (H5"D. the final result is: 



Ai-f = -eJ V-9(x) [/f (x) d^ f f {x)\ A»(x)d 4 x, (52) 

where the bilateral derivative acts as: g d h = {dh)g — (dg)h. 

The above expression is just the scattering amplitude for scalar field in the first order of 
the perturbation theory on de Sitter expanding universe. With the above calculations we 
establish the general rules of calculation, which can be also used in the language of Feynman 
graphs. The above reduction formalism and perturbation theory, allows us to deduce the 
correct definition for the scattering amplitudes. The results obtained in the last two sections 
will be applied to calculate the Coulomb scattering for a charged scalar field on de Sitter 
space-time. 

V. SCATTERING AMPLITUDE 

In the well establish picture from Minkowski quantum field theory, the scattering am- 
plitudes can be calculated without ambiguity, and the main results are the conservation of 
energy and momentum in scattering processes. We want to see how these results will be 
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affected by a curved background. This section will be entirely dedicated to the study of 
Coulomb scattering for a massive charged scalar field on de Sitter space. 

In our calculation we need the form of the Coulomb potential on de Sitter space-time 
which will depend on line element. Here the situation is immediate as we point out in [l6j, if 
one recalls the conformal invariance of the Maxwell equations and that the de Sitter metric 
is conformal with the Minkowski one. If we write S for the Coulomb field in the Minkowski 

\x\ 

space, then we find for the corresponding de Sitter potential: 

i 5 (x) = ^e- ui , (53) 

where the hated indices indicate that we refer to the components in the local Minkowski 
frames. We also observe that (1531) is just the expression from flat space with distances 
dilated/contracted by the factor e~ wt . 

Having the above elements we can proceed to develop the theory for the scattering am- 
plitude on de Sitter spacetime. Here the necessary requirements for develop the scattering 
theory is the global hyperbolicity of the space 9| and having a complete set of solutions 
of the free equation for incident field and scattered field (Born approximation), with the 
distinction between positive and negative frequencies. 

Let <pi{x) and <fif(x) be the freely propagating waves in the in and out sectors, and we 
assume that they are both plane wave of positive frequency. Our intention is to calculate the 
amplitude of Coulomb scattering for the external field (153]) and for initial and final states of 
the form: 

0/(z) = fp f (x), <Pi{x) = fa(x). (54) 

Replacing our quantities of interest ( f53l) and (J3J) in the definition of scattering amplitude 
( 1521) we observe that the dependence of t, x allows us to split the four dimensional integral 
into a spatial integral and a temporal one. Then is immediate that the spatial integral have 
exactly the form as in the Minkowski space [5j: 

e i(fi-Pf)S 4^ 

\x\ \ff - fi 

It is clear that the not so simple part is the temporal integral where is contained the influence 
of the gravitational field via the factor k. Our amplitude that contain the temporal integrals 
can be wreathed after we perform the bilateral derivative and use the recurrence relation 
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d 3 x n5 = — ^. (55) 



2J) as : 



-aZ 



8ir\p) -fi\ 



Pi 



r(2). 



r(i) 



(2), 



r(i) 



H^iPf^LMz) ~ Hlt\Pf*)HZi{PiZ 



dz 



+ P_l 
2 



c 



o 



r(2) 



r« 



(2) 



r(i)/ 



HHUiPfzW^P.z) - H% +1 (p f z)H£(p, 



dz 



(56) 



where a = e 2 and we pass to a new variable of integration : 

^.—ujt 



Z 



U 



(57) 



It is important to point out that in (156]) the integration limits for time variable corresponds 
to t — ±oo , we assume that the interaction extends into the past and future. 

The evaluation of the integrals is discussed in Appendix A, here we give the final result 
in terms of: delta Dirac function, hypergeometric Gauss functions, beta Euler functions and 
Heaviside step functions: 

— aZ 



A 



i-*S 



8ir\p} -fi 



;Bi 



where 



Bi 



^i±^-2i8(p f -p i ) + -e( Pi -p f ) 

y/PfPi Pt 



hi & 



Pf 
Pi 



.Pi 



+—9(p f -pi) 

Pf 



PfJ XPf 



The above formulas are our main result. In Eq. 



the notations are as follow: 



(58) 



(59) 



9k(x) 



S-y fl nJ >\ J t f\> 

e A 

sinh (irk) 

_| ± 

sinh 2 (7r/c) 

,ik 



2 F 1 (§, \ + ik; 1 + ik; x 2 ) 2F1 (§, § + ik; 1 + ik- x 2 ) 



B(\-ik;l+ik) 



B(-\-ik;l+ik) 



2 F 1 (§, I - ik- 1 - ik; X 2 ) 2^1 (5, | - ik; 1 - ifc; x 2 ) 



B(\+ik;l-ik) 



B(-\ + ik-l-ik) 



X 



sinh (irk) 

A. 



?J Pi (|, i + zfc; 1 + ifc; x 2 ) 2F1 (|, § + zfc; 1 + ifc; x 2 ) 



S(-i;l + ifc) 



flfi;l + ifc) 



.2' 



sinh (irk) 



2 F 1 (§, I - ik; 1 - ifc; x 2 ) 2-F1 (§, § - ifc; 1 - ifc; x 2 ) 



B{ 



l;i 



ifc) 



s(|; 1 



ifc) 



(60) 



16 



h k (x) 



e A 

sinh (ilk) 

+ X 

sinh 2 (7r/c) 

Ak 



2 F 1 (|,|-^fc;2-^fc;x 2 ) 2 
B(-\ + ik-2-ik) ' X 

2 F 1 (H+ik;2 + ik; X 2 ) 



■X 



2Fi{\,\-ik;-ik]x 2 ) 
B(-ik; \ + ik) 

2 F i (if + ^;^;x 2 ) 



X 



sinh (7rfc) 



sinh (yrfc) 



B(-\-ik;2 + ik) * B(\-ik]ik) 

2 F l Q, | + zfc; zfc; v 2 ) 2^1 (|, f + jfej 2 + *fc; x 2 ) 

B(|;ifc) S(-|;2 + ifc) 

2 iq (I, § - zfc; -ik- x 2 ) 2F1 (|, § - ifc; 2 - *fc; v 2 ) 



•X 



B(h-ik) 



2 ! 



£( 



-i-2 

2> Z 



life) 



•X 



(61) 



where we note y = — or — . 

The above formulas is our result and in what follows we will discuss some of their physical 
consequences. Before that we must state that the argument x i n relations (l60l) . (l6"Tl) must be 
considered in interval < x < 1 (note the argument in the 9 functions). This is the domain of 
convergence of the hypergeometric functions, because in the limit x ~^ 1 the hypergeometric 
functions diverge. This observation implies that in the terms that contain hypergeometric 
functions we must always consider that pf 7^ pi and from this we conclude that in contrast 
with the flat space case, here the law of conservation for modulus of momentum is lost. 

Let us see what happens when we make the Minkowski limit of (1551) . The contributions 
involving functions g^ (%) , hk (x) wn l vanish in the limit k — > 00 which is just the part of 
the amplitude responsible for momentum non-conserving case. Only the delta terms will be 
effective in ( l59l) . and we obtain: 

(Pf + Pi 



B n 



VPfPi 



-2i5(p f -pi). 



(62) 



We know that the corresponding factor which can be obtained studying the scattering of a 
charged scalar field on constant field Q/\x\ from the Minkowski case is: 

(E f + Ei 



B 



M 



%E~ 



-US{E f - E^ 



(63) 



where Ef, Ei are the Minkowski energies and E 2 = p 2 + m 2 . For obtain the ratio between 
de Sitter/Minkowski amplitudes it is useful to consider the ultra relativistic case (the case 
of large momenta) p> m,p ~ E ,and using (1531 . (162]) we obtain: 

26(E f -Ei) 2 y ' 

One expects that in the limit of a zero expansion of space uj — > 0, to obtain the Minkowski 
amplitude, since in this limit the metric and Klein-Gordon equation take the form from 
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Minkowski space. Our result may be a little surprising, but can be explained as follows. 
One finds after going back to the integration with respect to t that for t — > oo the integrand 
vanishes as e _a; ', while for t — ► — oo become divergent. As one see for t — > — oo the potential 
also become divergent. So the question is how come that we obtain a finite result, because 
a divergent behavior of the potential seems to be incompatible with a perturbational cal- 
culation from beginning. The answer might be as follows: as one goes back in time, with 
the divergent behavior of the potential, the contraction of the space will cause a similar 
divergent behavior in the incident frequencies as one see from (|3]), which lead finally to a 
finite result in our calculations. Using now the asymptotic formulas for Hankel functions 

A: 



when z takes large values 



H^ 2 \z) = (—} e ±i(*-"r/a-T/4) > s _+ oo, ( 65 ) 



and passing to conformal time t c = —z the contributions to the amplitude in (156]) will give 
( the contributions of all integrals from (56) is considered) : 

The integral in 1J661) produces the distributional term, integrating over an infinite interval. 
Observing that in f l66l) appears only the semi- infinite axis we obtain the 1/2 factor. To 
summarize, the fact that we don't recover the Minkowski amplitude is hidden in the expan- 
sion of space from infinite future because we don't have contribution from infinite future in 
( 166]) . One should not generally expect to recover Minkowski limit because the metric (JTJ is 
geodesically incomplete and is conformal only to a part of Minkowski space. 

The form of the amplitude (1551) . invites for some comments. First we observe that the 
result is a combination of terms with delta Dirac function and Heaviside step function, 
both having as arguments the difference between the final and initial (or initial and final) 
momentum. This means that the in wave which is a "monochromatic" plane wave with 
momentum p iy after scattering on our target is splitting in "monochromatic" plane waves 
which have the final momentum pf equal or different to the initial momentum p±. The 
probability of transition which is proportional to the squared amplitude is also a sum of 
terms that contain the two mentioned distributional terms. To summarize in de Sitter case 
the situation is different in comparison with the Minkowski case where we have a unique 
out wave, whose energy is equal with that of in wave and from here the conservation of 
momentum in the scattering process. 
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The result obtained in ([58]) for scattering amplitude lead to an important physical con- 
sequence which is the fact that the modulus of the momentum is not conserved in the 
scattering process. Let us make some comments about this result. Our conclusion is not so 
unexpected if one recalls that in Minkowski case the conservation of the momentum modulus 
is a direct consequence of the conservation of total energy which in turns is a consequence 
of the translational invariance with respect to time, which is lost in de Sitter space. 

VI. THE CROSS SECTION PROBLEM 

In the theory of quantum fields from Minkowski space the problem of calculating cross 
section is solved defining the probabilities of transition in unit of time. In our case the 
situation is complicated because our amplitude contain factors that broken the momentum 
conservation law. The terms from our amplitude that contain unit step functions don't 
allow us to define the probability of transition in unit of time as in flat space case. Only the 
delta Dirac terms 8{jpf — pi) allow us to define one probability of transition in unit of time 
-^ = -^7j4 (we attach the indices I to our probability to point out that the dependence 
of time is linear) which in our case give: 

dPi {aZf (p f +Pi) 2 x , ,d 3 p f 

nr = 32*»\p-,-M p fP , »'-% (67) 

For an explicit calculation of differential cross section we must evolve the incident flux, 
and for that we will use the formula for the scalar current of particles 

Jp = -i fp{x) d„ fp(x) , (68) 

from which we obtain for the spatial component: 

3(t) = -i [#(*) 5 fax)] = ^0^ 2) g e-') flg) (£ e-) . (69) 

The spatial components of the four vector density current of scalar particles is the density 
of incident flux. We see from (169]) that this is a quantity that depends of time which is 
unusual. 

According to the formalism of the scattering theory the amplitude obtained here has to 
be associated with the following picture: the in states is freely propagated from the moment 
t = back to the infinite past, then evaluated in presence of the interaction up to infinite 
future, and again freely propagated back to t — where is projected in out states. 
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For calculating the cross section in this complicated situation (incident flux dependent of 
time) we will make one approximation and follow the steps from [8|, where for calculating 
the incident flux we must know the state of the unperturbed system from the approximative 
moment of collision. Assuming that this moment is t ~ [8|, the incident flux will be: 

The differential cross section can be obtained easy because the scalar particles don't have 
spin: dai = tt-*- We obtain after using (IfJTj) . (1701) and performing the integration with 
respect to final momentum and use the fact that Pf — Pi — p , the differential cross section 
as follows: 



dai _ {aZfu r (2) /p\ (1) /p\i-i 



*VWZ\% 



(71) 



dn 64vrV sin 4 (f ) 

The Minkowski limit of the above cross section lead to a mathematical problem of great 
complexity which require a breakthrough in the theory of cylindrical functions, that is the 
approximation of Hankel functions when bought arguments have large values and one is 
imaginal and the other real. Our calculations in this direction show that this is a very 
difficult task, and lead to a series of difficulties, for that reason we can't give the Minkowski 
limit of (171]) . However one can observe that the discrepancy factor in the amplitude was 1/2 
and this lead us to the fact that the discrepancy factor when we take the Minkowski limit 
in (jnj) should be 1/4. 

In the case when the modulus of momentum is not conserved because the translational 
symmetry with respect to time is lost on de Sitter space-time, the probability of transition 
in unit of time can't be defined because the dependence of time is not linear. In this case the 
nonlinear part of amplitude at square has to be interpreted as the probability of transition 
(P n i ~ |A n ;| 2 ). Then the cross section in the non- conserving case is defined as the probability 
of transition P n i, from initial state into all states with momentum vector directed into dQ 
to the incident probability per unit area j (flux): 

da nl = ^ (72) 

3 

The probability of transition in the non-conserving case (dP n i = \A n i\ 2 j^p^) can be calculated 

using the properties of Heaviside functions and we obtain: 

{aZf 



dP 



ni 



647T 2 \p~f — Pi 



Pi \Pi / Pf 



■ P ' 
Pf k \Pf 



^- (73) 
(2tt)3' {76) 
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where the notations are as follows: 

X k {x) = \hk(x)\ 2 + Mx)\ 2 - 9k(x)hk(x) ' 9*k{x)K(xl (74) 

where x — Pf/Pi or Pi/Pf- Using (170]) and (1751) we obtain the cross section in the non- 
conserving case as: 



&nl 



a(oZ,! ^-sKr" 



32pi7r 3 \p f -p 



4 

1 1 



x / d 6 p f 



6{ Pl - Pf )- 2 X k ( V -l) + 9(p f - Pl )\x k (*l) 

pi \PiJ p f \PfJ 



(75) 



Pf 

As we can see from ( 1751) for obtaining the explicit expression of differential cross section an 
integration after pf is required. The dependence in functions gk{x)i hk{x) of final momentum 
is very complicated, and the integrals that must be solved contain products of hypergeometric 
functions and unit step functions. Solving the above mentioned integrals is not a easy task, 
their results are too complicated to be wreathed here , and we restrict to give the formulas 
in Appendix B. 

We see that our cross sections calculated in bought cases, depends on the form of the 
incident wave which is unusual. In Minkowski theory an essential role in the treatment when 
we deal with finitely extended wave, was played by the energy conservation factor S(Ef — Ei), 
which assure that the cross section will not depend on the form of incident wave. In our case, 
even in the momentum conserving case the cross section is dependent on the particular form 
of wave function of the incident particles. Thus seems that the cross sections calculated in 
the scattering processes on curved spacetimes will generally depend on the form of incident 
wave, which clearly is a consequence of the lost of translational invariance with respect to 
time. 

VII. CONCLUSION 

In this paper we construct the theory of scattering amplitude for scalar field on de Sitter 
space-time, using the solutions of the free Klein-Gordon equation which behave as plane 
waves with a given momentum. Our calculations shows that the methods used in flat space 
case for reduction of particles can be applied when we write the amplitudes of transitions 
in curved space-time. We obtain that every particle will be replaced after reduction with 
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expressions which depend on field operator 4>(x). The reduction method allows us to express 
the transition amplitudes as function of fundamental solutions of the free Klein-Gordon 
equation , and as function of vacuum expectation value of time ordered product of fields. 
Also we obtain that the vacuum expectation value of time ordered product of fields, together 
with normalization constants (\/z) define the Green functions of the interacting fields. With 
these elements we obtain the definition of the scattering amplitudes using the same method 
as in the flat space case. 

Further we investigated the Coulomb scattering amplitude for the charged scalar field in 
the expanding de Sitter space. We have complete ignored complications due to the ambiguity 
of the particle concept in a curved space-time. We considered the initial and final states of 
the field as exact solutions (with definite momentum) of the free Klein-Gordon equation in 
de Sitter space . Our calculations shows that the amplitude depends in an essentially way 
on the parameter k = y fi 2 — §■ The most important consequence of our result is the fact 
that the modulus of the momentum is not conserved in the scattering process on de Sitter 
space. In the limit of zero expansion of the space we don't recover the Minkowski amplitude, 
and we explain this as a consequence of the fact that our metric is geodesically incomplete. 
Another unexpected result is that the incident flux is a quantity that depends on time, and 
the cross section depends on the form of incident wave. All this unusual consequences are 
the result of the lost of translational invariance with respect to time on de Sitter space-time. 

For further investigations it will be interesting to study scattering processes in the second 
order of perturbation theory. However this will imply a series of technical difficulties related 
to the expression of Feynman propagators for Dirac field and scalar field on de Sitter space- 
time. Also another interesting subject related to the scattering processes is the construction 
of wave packets on de Sitter space, which could lead to a deeper understanding of the subject 
treated here. We hope that in a future paper we will present the above mentioned problems. 

APPENDIX A: INTEGRALS OF THE TYPE WEBER-SCHAFHEITLIN 



Here we give the main steps that 
formulas that we need are given in 



; help us to obtain the scattering amplitude fl58|) . The 
IB: 
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H^Xz) = '—^ z ' e ' TMt W z ) 



isin(7r/i) 



and the recurrence relation 



fl,0: 



2^^ = <-?(*)-<£>(*), (A2) 

which is used after we perform the bilateral derivation in our amplitude. Replacing the 
above formulas in ( l56i) we arrive at integrals of the type Weber-Schafheitlin 

a /T r + V +1 ) 



dzz s J ll (az)J l/ (bz) 



2 



^ u^ "<-l\jj,^t i)L y — ■— 



2»6^-«+ir(/x + i)r c / ~ /x + 8+1 ; 



_ f u + u — s + 1 u — u — s + 1 1 a 2 

x 2 F^ , ;„ + !;_ 

b > a > 0, i?e(s) > -1, i?e(/i + z/ - s + 1) > 0, (A3) 



and 



oo yn ^+y-s+l ' 



cfez s J lJL {az)Jy{bz) 



' pu — U+S + l ' 

- 2 ■ 



2 s a^~ s+1 r(z/ + i)r 

^ fa + v — s + 1 zz-u-s + 1 „ 6 2 

x 2-^i , — - ; v + 1; —z 

2 1 V 2 ' 2 ' ' a 2 

a > 6 > 0, i?e(s) > -1, i?e(/i + ^ - s + 1) > 0, (A4) 

we specify that the above two cases are fulfilled in our analysis. Our integrals demand 
s — — 1. We will evaluate the integrals in fl56l) directly with s — — 1 and with fi, v arguments 
following from (J2J). 

Let us take a look to the distributional contribution 5{pf — p^. The distributional term 
S(a — b) is generally present in the integrals ( lA3l) .f lA4l) for s = — 1 (the case of real arguments 



jjl — v is a familiar one). As we point out in 



16), here we can start with s = — 1 + e and 



integrate in ( 1A3I) . (1A4I) with respect to b over the interval a — e < b < a + e. Then the 
coefficient of 5(a — b) will be obtained then as the quantity which is left in the limit e — > 0. 
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APPENDIX B: INTEGRALS WITH HYPERGEOMETRIC FUNCTIONS 



For calculation of differential cross section (1751) .we solve integrals whose general form is: 

-I 2 



0(Pf - Pi 



2F1 I 6, c; d; -\ 
Pf 



pf k+n d Pf 



1 
Pi 



-ik—n/2 



P?r(d) 

2r 2 (6)r 2 (c) 



-r( c ) 



r(i - ik - §)r(6)r(i + 6 + «fc - §)r(i + c + «fc + f ) 



r(i + d + ik + %) 

Tt Tl Tt Tl 

X4F3 1 6, c, — 1 — ik , —d — ik ; d, —b — ik , — c — ik ; 1 

1 ' ' 2 2 2 2 



r(c)r(i + b + ifc + f)r(b)r(-| -ik- §)r 

r(i + d + ifc + f; 



n\r/ l+2c+2ik+n ~> 



"Pi A/ 2 

2 



(1 71 I 71 I 7Z 

b,c,---ik--,--d-ik--;d,--b-ik--, 

\ r _ . k _n \ r(rf)r(| + & + c + jfc + 1 )r(6)r(-| -b-ik- 1 )r(-6 + c) 

2~ C * 2' ; + r(-6 + d)r(i + d + 6 + ifc + §) 



,'l + 6 + rf + 2^ + n N 
xT ( | ^ 



1 / "I n 

4^3 ( 6, 1 + & - d, - + 26 + ik + 2 



1 77 i 77 1 77 

- + 6 + c + ifc + — ;l + & — c, - + 6 + ifc + -,- + 6 + cZ + ifc + — ;1) -r(-6 + c)r(6) 

^ Li La __ __ ^ 

r(d)r(-l -b-ik- f )r(l + 26 + iA; + f )r(l + 6 + c + ik + § ) 



r(-6 + d)r(l + 6 + d + iife + § ) 



n 



n 



n 



x 4 F 3 6, 1 + 6 - d, 1 + 26 + ik + -,1 + b + c + ik + -■ 1 + 6 - c, 2 + 6 + iA; + - 



n 



l + & + d + ifc + -;lH 



2' 2 ' 2 

r(c)r(| + 6 + c + %k + |)r(| + f + zfc + f )r(^)r( - 1 - 2c - 2 ^- n ) 



| / | 77, I 72 

xr(6- c)p iA / — g 4-^3 ( c, l + c-d,- + 6 + c + «£; + -,- + 2c + «£; + -; 1-6 + c, 
3 , n 1 , , n \ T(d)T(l + b + c + ik+ r j)T(-l-c-ik- r j) 

/ Th\ ( Tl 72 

xT \l + 2c + ik + -j aF 3 1 c, 1 + c-d, l + c + b + ik + -,l+2c + ik + -;l-b + c, 



n n 

2 + c + ik + -,l + c + d + ik + -;l 



(Bl) 



and integrals of the form , / °° 9(p i -pf) 2 F 1 \b, c; d; ^j 2 ^i (e, /; g\ ^j P/dp/; J °° 0(p/ 
Pi)2F\ (6, c;d;T-j 2-^1 ( e, /; #; t- J pjdpf , which is the same type as the above and have 
closer result. The above integrals are valid in the assumption n > 0. 
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